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We prove results on weak convergence for the alternating split Bregman algorithm in 
infinite dimensional Hilbert spaces. We also show convergence of an approximate split 
Bregman algorithm, where errors are allowed at each step of the computation. To be 
able to treat the infinite dimensional case, our proofs focus mostly on the dual problem. 
We rely on Svaiter's theorem on weak convergence of the Douglas-Rachford splitting 
algorithm and on the relation between the alternating split Bregman and Douglas- 
Rachford splitting algorithms discovered by Setzer. Our motivation for this study is 
to provide a convergent algorithm for weighted least gradient problems arising in the 
■ hybrid method of imaging electric conductivity from interior knowledge (obtainable by 

, MRI) of the magnitude of one current. 

ON 

1 Introduction 

Split Bregman and alternating split Bregman algorithms were proposed by Goldstein and 
Osher [8] for solving problems of the form 

x 

5h ' mm f(d) + g(u) subject to d — Du, 



where D £ f mxn is a linear transform acting from IR n to M m and /, g are convex func- 
tions. These algorithms have been shown to be very successful in various PDE based image 
restoration approaches and compressed sensing [HI [25]. In particular, the alternating split 
Bregman algorithm is very efficient for large scale / 1 -norm minimization and TV minimiza- 
tion problems problems [H |25] . The convergence of the split Bregman algorithm was proved 
in [8]. Later in [5] and [20| [2Tj the authors independently proved convergence of the alter- 
nating split Bregman algorithm in finite dimensional Hilbert spaces. In this paper we shall 
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prove weak convergence results for the alternating split Bregman algorithm in infinite di- 
mensional Hilbert spaces. Our motivation for this study is to provide a convergent algorithm 
for weighted least gradient problems arising in the hybrid method of imaging electric con- 
ductivity from interior knowledge (obtainable by MRI) of the magnitude of one current. Our 
proof relies on a recent result of Svaiter [22] about weak convergence of the Douglas- Rachford 
splitting algorithm. 

Let H\ and H2 be real Hilbert spaces and consider the minimization problem 

(P) mm{g(u) + f(Lu)}, 

where L : Hi — > H2 is a bounded linear operator and both functions g : Hi — >• RU {00} and 
/ : H2 — >■ K U {00} are proper, convex and lower semi-continuous. The problem (P) can be 
written as a constrained minimization problem 

min g(u) + f(d) subject to Lu = d, (1) 

u£Hi,d£H2 

which leads to an unconstrained problem: 

min g( u ) + f{d) + hLu-d\\ 2 . (2) 

u£H 1 ,deH2 I 

To solve the above problem, Goldstein and Osher [8] introduced the split Bregman method: 

d k+1 ) = argmin^^Mu) + f(d) + ± \\ b k + Lu - d \\*}, (3) 
b k+1 = b k + Lu k+1 - d k+l . 

Yin et al [21] (see also [23]) showed that the split Bregman algorithm can be viewed as an 
augmented Lagrangian algorithm [91 [161 EE]- Since the joint minimization problem ([3]) in 
both u and d could sometimes be hard to solve exactly, Goldstein and Osher [8] proposed 
the following algorithm for solving the problem (P). 

Alternating split Bregman algorithm: 

Initialize b° and d°. For k > 1: 



1. Find a minimizer u k of 



2 

on Hi. 

2. Find the minimizer d k of 



I k (u) := g(u) + - || b'- 1 + Lu- d k ~ l \\\, (4) 



I k (d):=f(d) + ^\\b k - 1 + Lu k -d\\l (5) 



on H2. 
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3. Let b k = b k ~ l + Lu k - d k . 

Cai, Osher, and Shen [5] proved that if the primal problem (P) has a unique solution then 
the sequence u k in the above algorithm will converge to the minimizer of (P). Independently, 
Setzer [20] showed that the alternating split Bregman algorithm coincides with the Douglas- 
Rachford splitting algorithm applied to the dual problem 



and proved the convergence of the alternating split Bregman algorithm in finite-dimensional 
Hilbert spaces (See [21], Proposition 1 and Theorem 4). 

We note that for general L, the functional I k {u) in the above algorithm may not have a 
minimizer in H\ and therefore the alternating split Bregman algorithm may not be well 
defined. We thus make the following definition, which will be used throughout this paper. 

Definition 1 We say that the alternating split Bregman algorithm is well defined if for all 
k > 1 the functionals I k {u) and I 2 {d) have minimizer s in H\ and H 2 , respectively. 

Note that if ( M ) is strictly convex and coercive in H 2 . Hence if / is weakly lower semi con- 
tinuous, then I k (d) will have a unique minimizer in H\. The following proposition provides 
a sufficient condition for existence of a minimizer of I k {u). 

Proposition 1.1 Let Hi and H 2 be two Hilbert spaces, g : Hi — > M U {00} be a proper, 
convex, and lower semi- continuous function. Assume L : Hi — )■ H 2 is a bounded linear 
operator. If L*L \ H\—t H\ is surjective, then for every c G H 2 the functional 



has a unique minimizer on Hi. 

We include a proof of Proposition 1.1 in Section 2 of this paper. We are now ready to state 
our main theorems. Below is a special case of a more general result (Theorem 2.6) that we 
will prove in Section 2. We let L* denote the Hermitian adjoint of L. Throughout the paper, 
we make the usual identification of a dual of a Hilbert space H with H itself. 

Theorem 1.1 Let Hi and H 2 be two Hilbert spaces (possibly infinite dimensional) and as- 
sume that both primal (P) and dual (D) problems have optimal solutions and that L*L : 
Hi — > Hi is surjective. Then the alternating split Bregman algorithm is well defined and the 
sequences {& fe }fcgv ond {d k }keN converge weakly to some b and d, respectively. Moreover Xb 
is a solution of the dual problem, {d k — Lu k+1 }teN converges strongly to zero, and 



(D) 



mm{g*(-L*b) + r(b)}, 




(6) 



CO 




(7) 



k=0 



Furthermore, there exists a unique u G Hi such that Lu = d and u is a solution of the primal 
problem (P). In particular {u k }k^N has at most one weak cluster point u = u. 
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Under additional assumptions on the functionals / and g we will also prove the following 
result. 

Theorem 1.2 Let Hi and H 2 be two Hilbert spaces (possibly infinite dimensional) . Suppose 
that alternating split Bregman is well defined, f is continuous, and both f and g are weakly 
lower semi continuous. Then the sequences {b k }k£N and {d k }k£N converge weakly to some 
b and d, respectively. Moreover Xb is a solution of the dual problem (D), {d k — Lu k+1 }k£N 
converges strongly to zero, and 

00 

^\\d k -Lu k+1 \\ 2 2 < 00. (8) 

k=0 

Furthermore 

lim f(Lu k ) + g(u k ) = min f(Lu) + g(u), (9) 

k— too udH\ 

every weak cluster point of {u k }keN is a solution of the problem (P), and L(u) = d for every 
cluster point u of {u k }keN- I n particular if (P) has a unique solution then the sequence 
{u k }k£N has at most one weak cluster point which is a solution of the problem (P). 

Comparing Theorem 11.21 with Theorem 11.11 (and with our more general result Theorem 12.61) 
one can see the effect of the operator L on the convergence behaviour of the alternating split 
Bregman algorithm. Indeed if L is injective then Theorem 12.61 guarantees that the sequence 
{u k }keN has at most one weak cluster point without assuming uniqueness of minimizers 
of the problem (P), continuity of /, or weak lower semi continuity of the functionals /, g. 
When L is injective and the problem (P) has more than one solution then, depending on the 
initial values of b° and d°, the alternating split Bregman algorithm may converge to different 
solutions of the primal problem. On the other hand when L is not injective and the primal 
problem (P) has a unique solution then Theorem 11.21 guarantees that {u k }k^N has at most 
one weak cluster point while Theorem 12.61 only says L~ 1 (d) contains a solution of the primal 
problem, which is a weaker conclusion. 

In many applications L*L is surjective and Theorem 11.11 guarantees the convergence of 
the alternating split Bregman algorithm. For instance in weighted least gradient problems 
Hi = #o(ft), H 2 = (L 2 (n)) n , and L = Vu. It is easy to check that L* : (L 2 (n)) n -> i? _1 (O) 
is surjective. Recently, our group studied the problem of recovering an isotropic conductivity 
from the interior measurement of the magnitude of one current density field EE2 EES] ■ We 
showed that the conductivity is uniquely determined by the magnitude of the current gener- 
ated by imposing a given boundary voltage. Moreover the corresponding voltage potential 
is the unique minimizer of the infinite-dimensional minimization problem 

M = argmin{ / |J||Vw| :veH l (Q), v\ 9n = f}, (10) 

where \J\ is the magnitude of the current density vector field generated by imposing the 
voltage / on the boundary of the connected bounded region Q C R n , n > 2. 
The results presented in this paper lead to a convergent split Bregman algorithm for com- 
puting the unique minimizer of the least gradient problem (jTt)]) . The details will be presented 
in a forthcoming paper [13], along with a number of successful numerical experiments for 
recovering the electric conductivity. 
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2 Convergence of the alternating split Bregman algo- 
rithm 



Recall that, by Fenchel duality [19] . the dual problem corresponding to the problem (P) can 
be written as 

(D) -min{g*(-L*b)+r(b)}. 

0G-H2 

Let v (P) and v (D) be the optimal values of the primal (P) and dual problem (D), respectively. 
Weak duality always holds, that is v(P) > v(D) [TJ [19]. To guarantee strong duality, i.e., 
the equality v(P) = v(D) together with existence of a solution to the dual problem, several 
regularity conditions are available in the literature (see [I], Chapter 7). In particular, if 

3x G dom(g) fl dom(/oL) such that / is continuous at Lx, 



or 



G int(Ldom(g) — dom(/)), 



then strong duality holds. In this paper we will always assume that both primal and dual 
problems have optimal solutions, and at least one of the above conditions is satisfied. Then 
by Rockafellar- Fenchel duality [19], if b is any solution of the dual problem, then the entire 
set of solutions of the primal problem is obtained as 

^(-L^nr^flfe). (11) 

The above representation of solutions of the primal problem is the key for our understanding 
of the alternating split Bregman algorithm. To explain this, we first prove the following 
simple lemma. 

Lemma 2.1 Let g : H\ — > KU {00} and assume L : Hi — )■ H 2 is a bounded linear operator. 
Assume one of the following conditions hold: 

1. There exists a point Lu where g* is continuous and finite. 

2. L* : H 2 — > Hi is surjective. 
Then 

d{g*o{-L*)){b) = -Ldg*{-L*b), (12) 
for all b G H 2 . In particular if L~ 1 {d) 7^ for some d G H 2 , then 

- d G d{g*o{-L*)){b) <=> L^d n dg*(-L*b) ^ 0. (13) 

Proof. If 1) holds, then (fl2j) follows from Proposition 5.7 in [7]. Now assume L* is surjective 
and let -L{u) = -d G d(g*o(-L*))(b). Then 

g\-L*{c)) - g*(-L*(b)) > (-Lu, c-b) = (u, -L*(c - &)), 

for all c G H 2 . Since L* is surjective, u G dg*(—L*b) and consequently L~ l d C dg*(—L*b). 
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Now assume u G d(—L*(b)). Then 

g*(-L*(c)) - g*(-L*(b)) > (it, -L*(c - b)) = {-La, c-b) = {-d, c-b), 
for all c G H 2 . Hence -d G d(g*o(-L*))(b). □ 

By above lemma, if one can find d G L(H{) such that 

dedf*{b) and -de%*o(-L*))(6), 

for some solution 6 of the dual problem, then 

L-\d) C dg*{-L*b) n (14) 

Consequently, by Rockafellar-Fenchel duality, every u G L -1 d will be a solution of the primal 
problem (P). One can find d G if 2 satisfying the above conditions by solving the inclusion 
problem 

Oed(g*o(-L*))(b) + df*(b). (15) 

Indeed if 6 G H 2 is a solution of the dual problem and d G df*(b), then — d G d(g*o(—L*))(b) 
and (j!4p follows from Lemma 12.11 This is summarized in the following lemma. 

Lemma 2.2 Let g : Hi — > R U {00} and / : H2 - > R U {00} fre two proper lower semi- 
continuous convex functions and assume that hypothesis of Lemma \2.1\ hold. Also suppose 
both primal (P) and dual (D) problems have optimal solutions and let b be an arbitrary 
solution of the dual problem, then every u G L _1 (9/*(6)) is a solution of the primal problem 
(P)- 

We thus focus on computing a solution of the problem ( |T5l) . This can be written in the form 
of an inclusion problem 

OeA(b) + B(b), (16) 

where A := d(g*o(—L*)) and B := df* are maximal monotone operators on H2. If both 
primal and dual problems have optimal solutions, then the above inclusion problem has at 
least one solution. Therefore, if we can find a solution b of the problem ( TLB"]) as well as 
d G B{b), then by Lemma [2T21 every u G L _1 (d) will be a solution of the primal problem (D). 
The Douglas-Rachford splitting method in Convex Analysis provides precisely such a pair 
(b,d). Following this route leads to the alternating split Bregman algorithm. Indeed it is 
shown by Setzer in |20j that the alternating split Bregman algorithm for the primal problem 
(P) coincides with Douglas-Rachford spliting algorithm applied to (TT6"j) (see Theorem 12.41) . 
We now explain this in more detail. Let H a real Hilbert space and A, B : H — > 2 H be 
two maximal monotone operators. For a set valued function P : H — )■ 2 H , let Jp to be its 
resolvent i.e., 

J P = (Id + P)- 1 . 

It is well known that sub-gradient of convex, proper, lower semi-continuous functions are 
maximal monotone [TS] and if P is maximal monotone then Jp is single valued. 
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Lions and Mercier [10] showed that for any general maximal monotone operators A, B and 
any initial element Xq the sequence Defined by the Douglas-Rachford recursion: 

x k+1 = (J A (2J B - Id) + Id- J B )x k , (17) 

converges weakly to some point x G H such that p = Jb{x) solves the inclusion problem 
(|T6|) . Much more recently, Svaiter |22j proved that the sequence p k = Js(x k ) also converges 
weakly to p. 

Theorem 2.3 (Svaiter [22J) Let H be a Hilbert space and A, B : H — > 2 H be maximal 
monotone operators and assume that a solution of [To]) exists. Then, for any initial elements 
xq and po and any A > 0, the sequences p k and x k generated by the following algorithm 

Xk+1 = JxA^Pk - X k ) +X k ~Pk 

Pk+1 = J\B{x k+1 ), (18) 
converges weakly to some x and p respectively. Furthermore, p = J\b{x) and p satisfies 

G A{p) + B{p). 

To apply Douglas-Rachford splitting algorithm one needs to evaluate the resolvents J\A{2p k — 
x k ) and J\B(%k+i) at each iteration. To evaluate the resolvents we are led to find minimizers 
of Ii (u) and l\ (d) in the alternating split Bregman algorithm. Indeed if we let 

A = d(g*o(-L*)), B = 8f*, x° = \{b° + d°), and p° = Xb°. 

Then the resolvents J\A(2p k — x k ) and J\B(xk+i) can be computed as follows 

J X A(2p k -x k ) = \(b k + Lu k+1 -d k ), 

and 

JxB(x k+ i) = \(b k + Lu k+l - d k+1 ), 

where u k+l and d k+l are minimizers of I k {u) and I^id), respectively (see [201 I2T] for a 
proof). The following theorem gives the precise relation between the sequences generated 
by the alternating split Bregman algorithm and those generated by the Douglas-Rachford 
splitting algorithm. 

Theorem 2.4 (Setzer |20j) The Alternating Split Bregman Algorithm coincides with Douglas- 
Rachford splitting algorithm applied to (D) with A := d(g*o(—L*)) and B := df* , where 

Xk = \(b k + d k ), p k = Xb k , k > 0. (19) 

The operator T := Ja(2Jb — Id) + Id — Jb is known to be firmly non-expansive, i.e., 
T — -Id + with R satisfying: 

|| Rx — Ry \\<\\ x — y \\ for all x,y G H. 

We will need the following lemma in our convergence proof. 
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Lemma 2.5 If T : H H is a firmly non-expansive operator and x k+ \ = T(x k ) with 
Xq G H , then 

II t II 2 _L II ll 2 ^ II "II 2 

|| x k+l ~ X || +|| X k +1 — X k || ^|| X}. — X || . 

Proof. Since T is firmly non-expansive, i? = IT — Id is a non-expansive operator. 
Hence 

|| ifc* - i?x || 2 = - || x k - x f +2 || Tx fc - Tx || 2 -2 || (Id - T)x k - (id - T)z || 2 . 
Therefore we have 

-(|| x k - x || 2 - || Rx k - Rx || 2 ) =|| x k - x || 2 - || x k+1 - x || 2 - || - x fe || 2 . 

Since R is non-expansive, the left hand side of the above inequality is non-negative, and this 
completes the proof. □ 

Now we are ready to prove our main theorem. 

Theorem 2.6 Let Hi and H 2 be two Hilbert spaces (possibly infinite dimensional) and as- 
sume that both the primal (P) and the dual (D) problems have optimal solutions. Suppose 
that ( fl^j) holds and the alternating split Bregman algorithm is well defined. Let {u k } k€ N, 
{d k } k£ N, and {b k } k ^N be the three sequences generated by the alternating Split Bregman algo- 
rithm. Then {d k } k ^N, and {b k } ke N converge weakly to some d, and b, respectively. Moreover 
Xb is a solution of the dual problem, the sequence {d k — Lu k+1 } ke ^ converges strongly to zero, 
and 

oo 

^||^-L M fc+1 || 2 < oo. (20) 

k=0 

Furthermore, L~ 1 (d) contains a solution u of the primal problem (P). In particular if L is 
infective, then {u k } k( zN has at most one weak cluster point u = u. 

Proof. The weak convergence of the sequences d k , and b k follows from Theorems 12.31 and 12.41 
To prove the estimate fl 3 . 3 H . let T = Ja{2Jb — Id) + Id — Jb- Since T is firmly non-expansive, 
by Lemma 12.51 we have 

II Xk+i - x || 2 + || x k+1 - x k \\ 2 <\\ x k -x || 2 , (21) 
where x is the weak limit of x k with T(x) = x. By the above inequality, we have 

oo 

^2 ii Xk + i ~ Xk n 2< 00 • ( 22 ) 

k=0 

Now observe that 

Xk - Xk _ x = \(b k+1 + d k+1 -b k - d k ) = \{Lu k+1 - d k ), 
and hence (13 .31) follows. 
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By Theorem 12.41 and Theorem 12. 3[ p = Xb is a minimizer of the dual problem and J\df* (X(d+ 
b)) = Xb. Therefore 

Xb + Xdf*{Xb) = X(d + b) <=> dedf*{Xb) & -d e d{g*o{-L*)){Xb). 

By Lemma [2.11 there exists u G Hi such that w G dg*(—L*(p)) and = d. Therefore 

uedg*(-L*(p))nL-\dr(p)). 

Since p is a minimizer of the dual problem, it follows from the Fenchel-Rockafellar duality 
theorem that u is a minimizer of the primal problem. If L is injective then the sequence 
u k has at most one weak cluster point u and necessarily u = u. The proof is now complete. □ 

Proof of Theorem II. 2t Since / is continuous by Theorem 4.1 in [7] the dual problem (D) 
has an optimal solution. Thus it follows from an argument similar to that of Theorem 12.61 
that the sequences b k and d k weakly converge to some b and d where Xb is a solution of the 
dual problem and (|S]) holds. In particular if u is a weak cluster point of u k then L(u) = d. 
To prove OH]), we can now use the argument of Cai, Osher, and Shen in the proof of Theorem 
3.2 in [5] Let u be a solution of the problem (P) and set d = Lu, p G df(d), and b — ?. 
Define 

u k = u k - u, d k = d k - d, b k = b k - b. 

Then, as in [5], 



^ * ii* - ii r + ii <% t - ii e +i in 

K K 

= Y,( d 9( uk+1 ) - dg{u),u k+1 - u) + Y,(9f(d k+1 ) - df(d),d k+l - d) 

k=0 k=0 

+ \ II Lu k ^ - rf^ 1 f + II - f) • (23) 

\fc=0 fc=0 / 

Now ([9]) follows from (123]) as in the proof of Theorem 3.2 in [5]. Finally since both / and g 
are weakly lower semi continuous, in view of OH]), every weak cluster point of {u k }keN is a 
solution of the primal (P). □ 



Remark 2.7 Notice that Theorem \2.cl\ is crucial for the proof of the convergence of the se- 
quences b k and d k . 



Proof of Proposition ll.lt First note that u is a minimizer of (Q if and only if 

G (1/X)dg(u) + L*Lu + L*c^ue {{l/X)dg + (-L*c). (24) 

Therefore to guarantee existence of a solution of ([6]) it is enough to prove that (1/X)dg + L*L 
is surjective. Let A := L*L. Since A is the subgradient of the convex lower semi-continuous 



9 



functional || Lu || 2 , it is a maximal monotone operator. We claim that A is 3*— monotone, 
i.e. 

V(x, x*) E Hi x Hi sup (x -y,y* - x*) < oo. (25) 

(y,V*)&Graph{A) 

Since A is surjective and (y,y*) G Graph(A), 

y* = A(y), and x* = A(xo) for some xq G Hi. 

Hence 

(x-y,y*-x*) = (L(x-y),L(y-x )) 

= -\\L(x - y)\\ 2 + (L(x - x ),L(y - x )) 

< -\\L(x -y)\\ 2 + \\L(x - x )\\\\L(y - x Q )\\ 

< i^z^U <00 . 

~ 4 

Thus (|25|) holds and A is 3*— monotone. It follows from Theorem 8 in [11] and Lemma 2.2 
in [3] that the operator {1/X)dg + L*L is surjective and therefore ([HD has a solution. 
Now notice that since L*L is surjective, L* is surjective and hence L is injective. Conse- 
quently the functional I(u) is strictly convex and has a unique minimizer. The proof is now 
complete. □ 

Remark 2.8 If A = L*L is surjective then one can show that A is also invertible. To see 
this note that if A is surjective then L* is surjective and hence L is injective. Now assume 
L*L(u) = 0. Then 

= {L*L{u),u) = \\Lu\\ 2 . 

Therefore u = 0. 

Theorem 11.11 now follows from Theorem 12.61 



3 Approximate alternating split Bregman algorithm 

In this section we show that the alternating split Bregman algorithm is stable with respect 
to possible errors at each step in the calculation of minimizers of If (u) and /|(c?). The proof 
relies on the following theorem about the Douglas- Rachford splitting algorithm. 

Theorem 3.1 (Svaiter [22J) Let A > 0, and let {ctk\k£N and {(3k}k<=N be sequences in a 
Hilbert space H . Suppose G ran(A + B), and ^fcgjvdl a k II + II 0k II) < oo. Take x G H 
and set 

Xk+l = %k + J~/A(2(J\BXk + Pn) ~ Xk) + OL k - (J\ B X k + 0k), k > 1. (26) 

Then x k and pk = J\BX k converge weakly to x G H and p G H , respectively and p = J\bx G 
(A + B)-\0). 
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The proof of the above theorem in infinite-dimensional Hilbert spaces is due to Svaiter [22] 
(see also jl]). 

Approximate alternating split Bregman algorithm: 

Initialize b° and dP. For k > 1: 

1. Find u k such that 

II Lu k -Lu k ex \\ H2 < a k , 

where u k ex is a minimizer of 

I*(u) = {g{u) + ~ || b k - l + Lu-d k - 1 

on 

2. Find such that 

II ^ — ^ez \\h 2 < Pk, 

where d k x is the minimizer of 

I k {d) = argmin^ 2 {/(rf) + ^ || ft*" 1 + Lu k - d || i>, 

on Hi- 

3. Let 6 fc = b k ~ x + Lw fc - d k . 

By Theorem 13.11 and an argument similar to that of Theorem 12.61 we can prove the following 
theorem about convergence of the sequences u k , d k , and b k produced by the above algorithm. 

Theorem 3.2 Let Hi and H2 be two Hilbert spaces (possibly infinite dimensional) and as- 
sume that both primal (P) and dual (D) problems have optimal solutions. Suppose that [IE) 
holds, the perturbed alternating split Bregman algorithm is well defined, and 

00 

^(ctfc + 0k) < 00. 

k=l 

Let {u k }k£N, {d k }keN, and {b h }k£N be the three sequences generated by the approximate 
alternating Split Bregman algorithm. Then {d h }k^N, and {b h }k^N converge weakly to some 
d, and b, respectively. Moreover Xb is a solution of the dual problem, the sequence {d h — 
Lu k+1 }keN converges strongly to zero, and 

00 

^||d fc -Lu fc+1 ||2<oo. 

fc=0 

Furthermore, L~ l (d) contains a solution u of the primal problem (P). In particular if L is 
injective, then {u k }k^N has at most one weak cluster point u = u. 
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If L*L is surjective, we have the following stronger result. 

Corollary 3.3 Let H\ and H2 be two Hilbert spaces (possibly infinite dimensional) and 
assume that both the primal (P) and the dual (D) problems have optimal solutions and that 
L*L : Hi — > Hi is surjective. Suppose that 

00 

^2(a k + (3 k ) < 00. 

k=l 

Let {u k }k£N, {d k }keN, and {b k }k&N be the three sequences generated by the approximate al- 
ternating Split Bregman algorithm, then the sequences {d k }keN, and {b k }keN converge weakly 
to some d, and b, respectively. Moreover Xb is a solution of the dual problem, the sequence 
{d k — Lu k+1 }keN converges strongly to zero, and 

00 

J2\\d k - Lu k+1 \\ 2 2 < 00. 

fc=0 

Furthermore there exists a unique u £ Hi such that Lit = d and u is a solution of the primal 
problem (P). In particular, {u k }k^N has at most one weak cluster point u = u. 
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